Full Paper 



Proc. of Int. Conf. on Control, Communication and Power Engineering 201 3 



Heuristic Optimization Technique for CHP- Wind 

Power Dispatch 

Vishal Chaudhary, Manjaree Pandit, Hari Mohan Dubey 
Department of Electrical Engineering 
Madhav Institute of Technology and Science 
Gwalior (M.P) India 

vishal.chaudhary30@yahoo.com; manjaree_p@hotmail.com; harimohandubeymits@gmail.com 



Abstract — This paper applies Differential Evolution (DE) for 
solving combined heat and power dispatch (CHP) problem with 
wind turbines. The CHP economic dispatch problem consist 
of two sub problems, namely, heat dispatch and power dispatch. 
The optimal dispatch utilization of multiple CHP units and 
wind turbines is a complicated problem due to the complex 
constraints such as feasible operating regions created by 
mutual dependence of heat and power. The uncertainty of 
wind speed and wind power is modeled using the Weibull 
probability distribution. Economic dispatch is carried out for 
cost optimization with varying wind power. The proposed 
method is able to handle all complex constraints very 
effectively without any violations. The developed optimization 
model is tested and validated on a study system from literature. 

Index Terms — Combined heat and power dispatch, Differential 
Evolution, Weibull probability density function, Wind power, 
Heuristic Optimization. 

I. Introduction 

Cogeneration or combined heat and power is simulta- 
neous production of multiple forms of energy in a single 
process. The Federal Energy Regulatory Commission (FERC) 
in the US defined cogeneration as "the combined production 
of electric power and useful thermal energy by sequential 
use of energy from one cause of fuel". The benefit of cogen- 
eration is achieved mainly from utilization of energy, which is 
normally lost in condenser during the separate power pro- 
duction. By capturing the excess heat, CHP reaches an effi- 
ciency of up to 89%, compared with 55% for the best usual 
plants. This means that a smaller amount fuel needs to be 
consumed to produce the same amount of useful energy. It 
also contributes to emission reduction compared to systems 
with separate electric power and heat production [ 1 ] .The pri- 
mary objective of economic dispatch in a conventional power 
plant is to find the optimal power outputs of generating units 
such that the total demand matches the generation with mini- 
mum fuel cost. In Ref. [2] a stochastic model for optimal dis- 
patch of CHP units along with wind power is carried out 
while treating harmful emissions as inequality constraints. 
In Ref. [3] a model of CHP dispatch is developed for a system 
consisting of both thermal generators and wind turbines. In 
this model, the chance of stochastic wind power is included 
in the constraint set and used the cross-product of power 
and heat terms, as well as two-dimensional polygons called 
feasible operation regions (FORs). Here Wind Power (WP) 
means the real electric power generated by wind turbines 
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and not the input wind power. Evolutionary optimization 
methods are very popular among researchers due to their 
fast convergence; population based effective search mecha- 
nism and simplicity. The major issue with evolutionary tech- 
niques is to create a balance between global search and local 
search. Issues such as premature convergence, setting up 
optimal tuning parameters and finding optimal population 
size appear prominently in recent publications [4-6]. One of 
the major advantage of wind energy is that, after the initial 
land and capital costs, there is basically no cost involved in 
the production of power from wind energy conversion 
systems(WECS).The primary feature that differentiates wind- 
powered from conventional generators in the economic dis- 
patch problem is the stochastic nature of wind speed [7] .Wind 
energy is an inexhaustible energy, which can be used to pro- 
duce electricity without the use of fuel and can save fuel cost 
for the power system. But on the further hand, wind energy 
is indiscriminate, which increases the insecurity during the 
operation of power system and presents new challenges to 
the economical operation of power systems. In Ref. [8] the 
economic dispatch problem for power systems which con- 
tain wind power generation plants. In Ref. [9] an extension of 
Simulated Annealing algorithm is developed for the solution 
of wind-thermal coordination scheduling problem. A signifi- 
cant part of installed thermal generation consists of com- 
bined heat and power (CHP) units, which enforce additional 
constraints on unit commitment and dispatch due to heat- 
demand-driven operation [10]. Wind is regarded as the most 
suitable renewable energy resource for bulk power genera- 
tion and the best alternative to the conventional energy re- 
sources mainly because of the increased capacities of wind 
turbine generators that are available [ 1 1 ] .The benefits of CHP 
and WP have been well known inside and outside the power 
engineering society[12].In Ref. [13] the problem of economic 
optimal operation of CHP system consisting of wind turbine, 
PV, fuel cell, waste heat boiler, gas boiler, battery, thermal and 
electric load is used. A plant growth simulation algorithm has 
been applied for optimal generation scheduling of power 
system having coal fired units and wind farm [ 14] . In place of 
traditional fossil fuel based power units the present trend is 
to incorporate non-conventional power generation through 
microgrid. Economic deployment of CHP based distributed 
energy resources considering cost as well as emission mini- 
mization is proposed using DE [ 15] . Electric heat pumps (EHP) 
can provide space heating in the domestic sector [ 1 6] . In Ref. 
[17] the effect of short duration wind generation variation on 
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optimal load dispatch is considered. DE is designed for dy- 
namic economic dispatch model including wind energy [18]. 

II. The Ed Model With Chp-Wind 

In this section, we establish an ED-CHP model 
characterized by the probability distribution. First, we define 
p to be the tolerance that the total demand cannot be 
satisfied. For example, if p =0.1, then up to 10% of the chance 
of insufficient supply could be tolerated. Therefore, a larger 
p implies more tolerance towardinadequate supply, and vice 
versa. Secondly, the model will take the transmission loss 
into account. For convenience of presentation, in the 
following, Total sum of total power demand (P d ) and total 
transmission losses (P ) is given by: 
P=P+P 

els d s 



represented by a small wind farm or a cluster of turbines in a 
large wind farm. Otherwise, more than one random variable 
should be included in the model. 
Constraint (4) can be rewritten as: 

Fw ~ X cpi ) = Fr ( w - p * ~ X "*) - * (7) 

Here F_, is the Cumulative distribution function of random 
variable W and p is the upper bound of probability that the 
sum of real power not greater than P ds . 

It should be mentioned that the value of p must satisfy 
the following relation: 



P (W=0)<p a <l. 



(8) 



= p i + Z Zt B ^ cp ^ cp i + X ***** + E w> 



(1) 



Where B is the parameters of transmission losses, cp. is 
the power generated by generator i, cp. is the power generated 
by generator j. 

With these notations, the ED-CHP model takes the 
following form: 

minimize 



This is because, if p <P (w=0), then the feasible region 
induced by (4) is *g- ^ m , Xj which implies that no WP 

would be pursued. Furthermore, p < 1 since p is a threshold 
of probability. 



V L=i 



(9) 



= 1 +■ expl- (— )*] 



" = ^fe ^ + feffl + % + + *il<*f + MM © - esp{-^- [v m » T + (v, - vJiP^ - £ eft)]*} 



subject to 



(3) 



Here v , v. ,v are the rated, cut-in, and cut-out wind 

r in out 

speeds, c is the Scale factor of the Weibull distribution, and 
k is the Shape factor of the Weibull distribution 

The above inequality can be easily converted into the 
following expression: 



Here q ,q. If q are the Power cost coefficients of generator 
i and P is the probability, the coefficients h. Q , h ,h are the 
heat cost coefficients of generator i and z is Cost index in the 
ED-CHP model. 

Heat generated by the generator i is given by 



X^ 

L =1 



^ ^ = Hi 



V — V; 



■ {— ln[ 1 + exp 



Where H is the total heat demand 



< bf„, 



(4) 



(5) 



(6) 



(10) 



Here i= 1 , 2, , N represents number of generating units 

and A is the 2N x 2N coefficient matrix of feasible operation 
region (FOR), b f is the 2N xl coefficient vector of FOR, cp is 
the N xl vector of power generation, ch is the N xl vector of 
heat output of generating units,, W represents all wind power 
(WP) to be dispatched. Since the total WP is characterized 
by a single random variable here, it implies that all wind 
turbines are located in a coherent geographic area, 
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Here w Rated power of the WP generating unit. 
Where Intermediary parameter (h) is given by 

C / V^ t \ i Vj„ 

P «V - *in V c* / v r - v-„ 

Therefore from (10) constraint (4) can be written as 



(11) 



(12) 
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III. Stochastic Wind Power Model 

Weibull distribution [19] is used for modeling the highly 
nonlinear relation between the wind power W and wind speed 
V (m/sec). It varies randomly with time and can be expressed 
in terms of the cumulative probability distribution function 

as: 



exp [-(-}*] 
c 



(v>0) 



(13) 



Where c and k are the scale factor and shape factor, 
respectively. Correspondingly, the probability density 
function of V is: 



k v . , v , 

fM = -(-y^exp [-(-)"] 
c c c 



(14) 



The relation between the input wind power and the output 
electric power depends on many factors, like the efficiencies 
of generator, wind rotor, gearbox, and inverter, depending on 
what kind of turbine is being employed. The generic relation 
between wind power and wind speed can be given as [20]: 

W= 







(V ~ *'m>> 



V|"v 



(V < tfr or V > v ear ) 
( v r < V < v 0ilC ) 



(15) 



v < v r ) 



Where v , v. , v : Rated, cut-in, and cut-out wind 

r in out 

speeds; W is the power generated by the wind generating 
unit; Wr is the rated power of the wind generating unit; c is 
the scale factor and k is the shape factor of the Weibull 
distribution. 

IV. Differential Evolution 

Heuristic optimization techniques based on the principle 
of evolution have assumed marked popularity over the con- 
ventional methods of optimization, particularly for practical, 
constrained engineering problems. These methods are highly 
advantageous due to their powerful convergence character- 
istics and ability to generate feasible solutions for all kinds 
of problems. Differential evolution (DE) proposed by Storn 
and Price is a population based technique which uses ran- 
dom operators and real number. The DE algorithm employs a 
special differential operator to generate new solutions from 
the parent population in place of the classical crossover or 
mutation operators used in the well established genetic algo- 
rithm (GA) technique. The difference between two randomly 
selected parameter vectors is used to improve the popula- 
tion. The DE has three powerful stages named as mutation, 
crossover and selection. 

A. Initialization phase 

An initial population of feasible solutions is generated 
randomly to represent the search domain using the following 
expression. 

X y =X-+rWx(X- -Xf )forj=l,2....D (16) 
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Where D is the dimension of the population, rand is an 
operator which generates a random number between 0- 1 and 
X"™ and X max are the minimum and maximum limits of the i* 
variable. The j" 1 member of the D dimensional population 
vector of size R can be represented as 



] fori=l,2 R (17) 



B. Mutation 

The initial target population is employed to produce a 
mutant/donor vector as shown below. 

Z,(f + 1)= X irl (t) + F[X ir2 (t)- X ir3 (t)] (18) 

The mutant vector Z is formed using two random 
population members other than X . The mutation rate F 
controls the differential variation that is instrumental in 
creating a mutant population to avoid search stagnation. It is 
selected between 0-2. Selection of F has a significant effect 
on the performance of DE algorithm. A smaller mutation rate 
makes the algorithm slow and longer time will be taken for 
convergence. A large F provides high level of exploration, 
but may cause the algorithm to skip optimal solution. Actually, 
the value of F should be small enough to enable the algorithm 
to explore tight valleys and large enough to allow global 
exploration in order to maintain diversity. The optimal value 
of F can be found by trial and error. 

The manner, in which the mutant vector is formed, 
depends on which DE scheme is employed. In this paper the 
mutation strategy known as DE/rand-to-best/1 as given below 
is found to perform the best. 

2#+l) =^(0+iT% e iO-^(0]+^, rl (0-^ >2 (0] (19) 

Here r l and r, are random integers (and not the same as the 
integer i ) generated in the range 1 -R 

C. Crossover Operation 

After the mutation operation is completed as explained 
above, the crossover operation is performed. In this operation 
some variables of the target vector (X.) are replaced by the 
mutant vector (Z.) to form the trial vector (U.) using the 
following logic. 



ujt+iy- 



[Z s (t + 1), . . ..if(Nun(j) < CR)orij = randintii)) 
]X S (f), if(Nunij) > CR)otij * randint(i)) (20) 

Where Nam (j) is a random number generated for the j th 
dimension and CR is a crossover rate selected in the range 
[0, 1]. 

D. Selection 

Selection is an important step and plays a major role in 
producing better solution/offspring. To decide whether or 
not the trial vector should be a member of the population of 
the subsequent generation, it is compared with the 
corresponding target vector. The selection process for an 
objective function under minimization is carried out by 
comparing the cost of each trial vector U.(t+1) with that of 
its parent target vector x.(t) as follows: 
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Viit + 1). ,.. if (u(t + 1)} < nxXt)) 

X,_{t). Btksrwise (21) 



V. Numerical Results And Analysis 



The CHP-wind model described by (2)-(6) a DE algorithm 
was developed in MATLAB environment. The constraints 
were handled by using penalty function approach. The 
algorithm was tested on a standard test function from 
literature [3] and the results are found to be better than the 
reported results. 

A. Description of the Test System 

The test system has two thermal /heat generating units 
and one wind farm. The cost calculation for unit one and two 
using coefficients for power and heat is given as C and C,. 



q = 0.1035 *cp, 2 +34.501 l*cp, +2650.03 12 + 0.0249&7j, 2 
+ 2.203 1*c/j, + 0.05 1 * cp, * c\ 

C 2 = 0.0435 * cp\ + 36.00 12* cp 2 + 1 250.043 1 + 0.0270chj 
+ 0.6002 *ch 1 + 0.0409 * cp 2 * ch^ 

The object is to minimize (Cj+C,) subject to the unit 
operating limit constraints given by 

6< cp, < 40; 0; ch, e 55 and 



(22) 



(23) 



10£c P] 
£ ch. 



60; 
; 30 



The FOR constraints of CHP units are given by 
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These constraints given by (5) get converted into straight 
line boundaries and the enclosed region represents the FOR 
of the two units, given as 

Jcpi + 14cAi < 840 
Vispi -Qch ± < 680 
-7rpi - 20c^ < -240 
cp 2 +■ 3ch 2 < 120 
Scp 2 + -2cb^ < 90 
Power Demand and Heat Demand are 80 MW and 70 
MWth. The cut-in, cut-out and rated wind speeds are 5, 45 
and 15 m/s respectively. The rated wind power is 20MW. 
Transmission losses are neglected in this study. 

B. Effect of probability of wind power 

The optimal dispatch of CHP units, i.e. power and heat is 
computed for different values of wind probability. The results 
are plotted in Table I. The results show that the cost of thermal 
units reduces with increasing wind probability. The Fig. 1 
plots the relationship between probability and wind power 
for different values of shape factor. It can be seen that the 
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wind power and wind reliability have a direct relation. Fig. 3 
shows the stable convergence characteristics of the De 
algorithm for CHP-wind problem for different population sizes. 
The mutation and crossover rate in DE were taken as 0.8 and 
0.7 respectively because best results were obtained at these 
values. 

Table I. Optimal Solutions (K=1.5) (PD=80MW, HD=70 Mwth ) 



n. 


0.25 


0.30 


0.35 


0.40 


0.45 


0.50 


0.55 


cpi 


47.1 


45.1 


43.1 


41.0 


38.9 


36.7 


34.4 


chi 


40.0 


40.0 


40.0 


40.0 


40.0 


40.0 


40.0 


q> 


30.0 


30.0 


30.0 


30.0 


30.0 


30.0 


30.0 


ch 2 


30.0 


30.0 


30.0 


30.0 


30.0 


30.0 


30.0 


Wind 


2.8 


4.9 


6.9 


8.9 


11.0 


13.2 


15.5 


Cost($/h) 


7179.3 


7086.4 


6994.2 


6901.2 


6806.5 


6709.1 


6607.7 


Power 
violation 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


Heat 
violation 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


Si . 1 1 1 1 1— 1 • 1 




Figure 1 . Variation of wind power with probability factor 
C. Effect of shape factor 

The shape factor k of Weibull distribution directly affects 
WP availability and hence cost of a wind-CHP system. Table 
I, II and III give the complete dispatch results of the test 
system for different values of shape factor. With increase in 
the value of shape factor the wind availability increases and 
therefore the cost of conventional power generation goes 
down. Fig. 2 shows the relation of generating cost and 
probability for different shape factors. 

Table II. Optimal Solutions (K=2.0) (PD=80MW, HD=70 Mwth ) 



P., 


0.25 


0.30 


0.35 


0.40 


0.45 


0.50 


0.55 


cpi 


43.9 


42.0 


40.3 


38.5 


36.8 


35.0 


33.1 


ch, 


40.0 


40.0 


40.0 


40.0 


40.0 


40.0 


40.0 


q> 


30.0 


30.0 


30.0 


30.0 


30.0 


30.0 


30.0 




30.0 


30.0 


30.0 


30.0 


30.0 


30.0 


30.0 


Wfad 


6.0 


7.9 


9.6 


11.4 


13.1 


14.9 


16.8 


Q)st(Vh) 


7030.5 


6947.6 


6867.6 


6789.3 


6711.5 


6633.2 


6553.4 


Fbwer 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


violation 
















Feat 


0.0 


0.0 


0.00 


0.0 


0.0 


0.0 


0.0 


violation 
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Table III. Optimal Solutions (K=2.5) (PD=80MW, HD=70 Mwth ) 



x 10 v 



ft 


0.25 


0.30 


0.35 


0.40 


0.45 


0.50 


0.55 


cpi 


41.7 


40.1 


38.5 


37.0 


35.5 


34.0 


32.5 


chi 


40.0 


40.0 


40.0 


40.0 


40.0 


40.0 


40.0 


cp> 


30.0 


30.0 


30.0 


30.0 


30.0 


30.0 


30.0 


Ah 


30.0 


30.0 


30.0 


30.0 


30.0 


30.0 


30.0 


Wind 


8.2 


9.8 


11.4 


12.9 


14.4 


15.9 


17.4 


Cost($/h) 


6933.4 


6859.1 


6790.1 


6723.1 


6657.4 


6592.3 


6526.8 


Power 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


violation 
















Heat 


0.0 


0.0 


0.00 


0.0 


0.0 


0.0 


0.0 


violation 

















D. Constraint handling and convergence behavior 

The idea behind the CHP-wind optimization is to utilize 
available wind generation and minimize the following 

function. The constants a and f3 are the power and heat 

penalty coefficients respectively. The second term imposes 
a penalty on the individual in terms of increased cost, if power 
balance constraint is not satisfied and the third term imposes 
a penalty, if heat balance constraint is not satisfied. 



Q+C 7 



Y^c Pj +WP-P Ds 



7=1 



Yfhj-H D 



7=1 



(24) 




Figure 2. Variation of cost with wind probability factor 
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Figure 3. Convergence characteristics of DE method for wind- 
CHP plant 

E. Comparision with Published results. 

The results of the DE based proposed method for optimal 
CHP-Wind dispatch is compared with results from literature 
[3]. The results obtained by the proposed method are found 
to be very close to [3]. Table IV presents the comparison of 
the two methods. The effect of Shape factor of the weibull 
distribution on cost of operation is studied. It is observe that 
as the value of k increases, the optimal cost decreases. For 
k=1.5 the cost was 7179 which reduces to 6933 when k=2.5. 
Similarly when the probability of wind power is increased the 
cost is found to reduce significantly 

Conclusions 

The proposed method is capable of computing feasible 
results and the proposed algorithm has a stable convergence 
behavior. The results are validated and compared with avail- 
able results in literature. The paper explores the relationship 
between the stochastic wind power and cost of a CHP-wind 
system. The selection of probability and shape factor for 
wind availability prediction play an important role in design- 
ing and planning a combined wind-CHP generating plant. 



Table VI. Comparison of cost 





k=1.5 


k=2.0 


k=2.5 


Pa 


0.25 


0.3 


0.35 


0.45 


0.25 


0.3 


0.35 


0.45 


0.25 


0.3 


0.35 


0.45 


Prop 
osed 
Meth 
od 


7179.3 


7086.4 


6994.2 


6806.5 


7030.5 


6947.6 


6867.6 


6711.5 


6933.4 


6859.1 


6790.1 


6657.4 


[3] 


7179 


7087 


6994 


6807 


7031 


6948 


6868 


6712 


6933 


6860 


6790 


6658 
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